(AD-AO4 lo G22) 
CETA 77.1 


A Simple Computer Model 
for Evaluating Coastal Inlet Hydraulics 


es (7 Wop u 


William N. Seelig 


Ze 
Sa 


DO LUMEN MF 
COASTAL ENGINEERING n=" Ton 7 
TECHNICAL AID NO. 77-1 ie 


JULY 1977 


ebro ed fer public release; 
distribution unlimited. 


U.S. ARMY, CORPS OF ENGINEERS 
a COASTAL ENGINEERING 


hs RESEARCH CENTER 
- Kingman Building 


Fort Belvoir, Va. 22060 


vaterial shall give appropriate 
ng Research Center. 


Contents of this report are not to be used for advertisin, 


authorized docu 


UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


REPORT DOCUMENTATION PAGE 


1. REPORT NUMBER 2. GOVT ACCESSION NO.) 3. RECIPIENT'S CATALOG NUMBER 


CETA 77-1 


4. TITLE (and Subtitle) 5. TYPE OF REPORT & PERIOD COVERED 
A SIMPLE COMPUTER MODEL FOR EVALUATING COASTAL Constal Bncimaeetine 
INLET HYDRAULICS Technical Aid ‘ 


6. PERFORMING ORG. REPORT NUMBER 


7. AUTHOR(s) 8. CONTRACT OR GRANT NUMBER(s) 


William N. Seelig 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 10. PROGRAM ELEMENT, PROJECT, TASK 
Department of the Army 


Coastal Engineering Research Center (CERRE-CS) F31019 
Kingman Building, Fort Belvoir, Virginia 22060 


AREA & WORK UNIT NUMBERS 


- CONTROLLING OFFICE NAME AND ADDRESS 12. REPORT DATE 

Coastal Engineering Research Center 13. NUMBER OF PAGES 

Kingman Building, Fort Belvoir, Virginia 22060 
- MONITORING AGENCY NAME & ADDRESS(i/f different from Controlling Office) 15. SECURITY CLASS. (of this report) 


UNCLASSIFIED 


1Sa. DECL ASSIFICATION/ DOWNGRADING 
SCHEDULE 


Approved for public release, distribution unlimited. 


DISTRIBUTION STATEMENT (of this Report) 


DISTRIBUTION STATEMENT (of the abstract entered in Block 20, if different from Report) 


- SUPPLEMENTARY NOTES 


KEY WORDS (Continue on reverse side if necessary and identify by block number) 


Coastal inlet hydraulics Computer model Computer program 


ABSTRACT (Continue am reverse sidn if necessary and identify by block number) 

A computer program for the prediction of coastal inlet velocities, 
discharge, and bay level fluctuations is presented. Two examples are given 
to demonstrate the numerical model. The computer documentation is included 
as an appendix, and the card deck may be obtained at CERC. 


DD Nong 1473 ~—s EprTion OF 1 NOV 6515S OBSOLETE UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


a ee ope alana sonseioney Sate eT ee 
Lay ee ind) AAS, 
reve 
hs ee Kat 
, } 


APTS a MER FUR. 

let ea re 

ae ee emotes dies mereLAy spt, sony - RE ny a: AEST AS i i 
a ; : 


Hipte 1 w & rhe oh) Gat 


sity ee a nt Bas eM } Tie 
a eae aie ek meee Fs Atl RR ) ‘a i : bi 
' a eid THe ee Saeuiene ry Ae oe Ae) a 

a 3 ee 


io ‘ Pie! a a 4 q =, 
sg Sa FARRER eee ss) (alarms a 1 i?) 


ae ee 
| ; Vaays | eee Seca yey pees 


aa), 8d eapOw Tat: Hain 
ee wee ue Sanaa 
F 4 wit - ats Tee hie ea 


oe eh APART MAR: seve Tenses! yici te 
Bel i tte LR Wd as oe oo ava : 


! 
» - - —_ dheien a 
H y 
' pik os 1h ow 
\ ns e amet ee Ae Mas Ran igs ab Ba No thet. 
i ian. qe) lietiw le Ph me, heel an ere tae An tain ait \w) 
i : 
} SIsSies as tins Me i 
- #— wAbee —— ww ee ee ee nc et lo ny 
‘ 
Sd ‘ 
ey 
) 
7 ; ies p 4 7 
-. “ — ferent = _— . aareayaris i haem wren 
; 7 a met yore Tarinase SS Wels: ob hie) me 
4 ‘, f 
, i ‘ hee Di : Pes Fy 
j way bebean MG) DR Raye th 3 
i\y ns bd 
% Shu 
7 j j ’ ‘a 


4 : g ‘ 
; paul : : aes Carcenyieried BLL alee ee, a 
[ . a oe Ye vol GM Shr cha: 
o “a isTni tagged Wik be Ae) phy’ tats i 
| a) Ts = ane Owl , yt PyiA ert PR wer uh 1 (iw A k 
t 


i ; a 1 “ik keV hime! Lins ni Me | bod 6 ae Dal 
i 


y fa Dr cea ; 1 aie 4 we Chand ee hi ’ 
7 a re , ae 
, ‘ ine oy 
| 
i 
| 
he: ——— _——a 


a VR PT en nie Tey 


PREFACE 


This report describes a method for estimating inlet velocities, 
discharge, and bay levels based on the numerical model of Seelig, Harris, 
and Herchenroder (in preparation, 1977). This method for predicting 
inlet hydraulics is not discussed in the Shore Protection Manual (SPM) 
(U.S. Army, Corps of Engineers, Coastal Engineering Research Center, 
1975). The work was carried out under the General Investigation of 
Tidal Inlets (GITI) of the U.S. Army Coastal Engineering Research 
Center (CERC). 


The report was prepared by William N. Seelig, Research Hydraulic 
Engineer, under the general supervision of Dr. R.M. Sorensen, Chief, 
Coastal Structures Branch. 


Comments on this publication are invited. 


Approved for publication in accordance with Public Law 166, 79th 
Congress, approved 31 July 1945, as supplemented by Public Law 172, 
88th Congress, approved 7 November 1963. 


OHN H. COUSINS 
Colonel, Corps of Engineers 
Commander and Director 


CONTENTS 


CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (SI) . 


SYMBOLS AND DEFINITIONS ..... 


I TNARODUCIUONS 66 6 66 6 56 6 6 6 


Teli PREDICTING INLET HYDRAULICS ... 


1. Systems Modeled with Computer Program. 
2. Procedures for Use of Computer Program 


IEIEIL EXAMPLES OF COMPUTER PROGRAM PREDICTION 
1. Cabin Point Creek, Virginia. ... 
2. Pentwater Inlet, Michigan. .... 
IV SWHMINR 6-6 6.6 6 6 0.0 61610 6 0 05 6 


TABLE 
Predicted Cabin Point Creek hydraulics. 


FIGURES 
IMILGESIDEY SWSEIN6 6 o o 6 6 6.06 070 0 


Calon iPelme Cisse, Walteeatmale: 6 59 5 6 0 0 
Cabin Point Creek cross-section... . 
Cabin Point Creek sea and bay levels. . 


Pemeanenecre Inulee 5 Whole 6 6 6 6 6 6 0 


° 


APPENDIX COMPUTER PROGRAM DOCUMENTATION (INLET). 


° 


Pentwater Inlet model prediction of monochromatic 


Pentwater Inlet model calibration... 


9 


forcing 


Predicted Pentwater Inlet velocities, discharge, and bay 
levels, and relative magnitude of terms in the equation 


OLAMOELOM. feornaemueuiccmen tes Felinccltcslinoiaonirs 


Page 


20 


CONVERSION FACTORS, U.S. CUSTOMARY TO METRIC (ST) 
UNITS OF MEASUREMENT 


U.S. customary units of measurement used in this report can be converted 


to metric (SI) units as follows: 


Multiply by To obtain 
nn ee eee —______ 
inches 25.4 millimeters 

2.54 centimeters 
Square inches 6.452 square centimeters 
cubic inches 16. 39 cubic centimeters 
feet 30. 39 centimeters 
0. 3048 meters 
Square feet 0.0929 square meters 
cubic feet 0.0283 cubic meters 
yards 0.9144 meters 
Square yards 0.836 square meters 
cubic yards 0.7646 cubic meters 
miles 1.6093 kilometers 
square miles 259.0 hectares 
knots 1.8532 kilometers per hour 
acres 0.4047 hectares 
foot-pounds 1.3558 newton meters 
millibars LOI, x 10-2 kilograms per square centimeter 
ounces 28.35 grams 
pounds 453.6 grams 
0.4536 kilograms 
ton, long 10160 metric tons 
ton, short 0.9072 metric tons 
degrees (angle) 0.1745 radians 
Fahrenheit degrees 5/9 Celsius degrees or Kelvins! 


eS eS OOOO OOO eeeeeeeeeeeeeees=<=~$~$—0M909M@M9090»S—SSsSSS\—<woooomo 
ITo obtain Celsius (C) temperature readings from Fahrenheit (F) readings, 


use formula: C = (5/9) (F -32). 


To obtain Kelvin (K) readings, use formula: K = (5/9) (F -32) + 273.15. 


SYMBOLS AND DEFINITIONS 
bay surface area (square feet) 
bay surface area at datum (square feet) 
coefficients to evaluate Manning's n (dimensionless) 
depth of bay (feet) 
maximum water depth in inlet (feet) 
Stillwater depth (feet) 
acceleration of gravity (32.2 feet per second squared) 
water level in bay (feet) 
water level in sea (feet) 
length of bay (feet) 
length of inlet (feet) 
forcing wave period (seconds) 
time step used in model (seconds) 


bay surface area variation parameter (dimensionless) 


A SIMPLE COMPUTER MODEL FOR EVALUATING COASTAL INLET HYDRAULICS 


by 
Willtam WN. Seeltg 


I. INTRODUCTION 


This report describes a method for estimating coastal inlet veloc- 
ities, discharge, and bay levels using the simple numerical model of 
Seelig, Harris, and Herchenroder (in preparation, 1977)!. The model 
can be used for sea level fluctuations caused by astronomical tides, 
storm surges, seiches, or tsunamis. A digital computer program is used 
because of the large number of computations. A run on a CDC 6600 com- 
puter generally costs less than $5 for a tidal cycle. 


II. PREDICTING INLET HYDRAULICS 


1. Systems Modeled with Computer Program. 


An inlet-bay system consists of a "sea" (e.g., ocean or lake) con- 
nected to a "bay'' by one or more inlets (Fig. 1). The computer model 
will predict bay levels, inlet velocities, and discharge as a function 
of time given the geometry of the system and the water level fluctua- 
tions in the sea. It is assumed that the sea is much larger than the 
inlet and bay and that the bay is large compared to the inlet. 


The model is designed for systems where the bay water level rises 
and falls uniformly throughout the bay. This occurs when the wavelength 
in the bay is much longer than the longest axis of the bay: 


Tp Y8dpay >> “bay > (1) 
where 
Ty = forcing wave period 
g = acceleration of gravity 
db ay = depth of bay 


Lp ay = length of bay 


2. Procedures for Use of Computer Program. 


Step 1. Evaluate the inlet geometry by using maps, charts, hydro- 
graphic surveys, and dredging records to determine the depth of water 
throughout the inlet. The side slope of the inlet at mean water level 


ISEELIG, W.N., HARRIS, D.L., and HERCHENRODER, B.E., '"'A Spatially 
Integrated Numerical Model of Inlet Hydraulics," GITI Report 14, U.S. 
Army Corps of Engineers, Coastal Engineering Research Center, Fort 
Belvoir, Va., and U.S. Army Engineer Waterways Experiment Station, 
Vicksburg, Miss. (in preparation, 1977). 
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should also be measured. Whenever possible, obtain this information for 
the time of interest because inlets frequently change shape, especially 
during major storms. 


Step 2. Construct a flow net (series of cross sections and channels) 
for the inlet to represent the model grid (Fig. 1). The flow net and 
inlet discharge are used to determine bottom friction throughout the in- 
let. The flow net is drawn by approximating the average path (channel) 
that water follows during ebb flow and floodflow. Channel boundaries 
are drawn along these paths for up to seven channels. A simple inlet 
with constant depth and width may be modeled with one or two channels. 
Complex inlets require approximately three to seven channels. Channels 
should have the smallest spacing in deep parts of the inlet where flow 
will be highest. Up to eight cross sections should then be drawn per- 
pendicular to the channels. The first cross section in the sea and the 
last cross section in the bay should have cross-sectional areas 10 times 
larger than the minimum cross-sectional area. Cross sections should be 
drawn with the narrowest spacing near the minimum cross-sectional area 
section where friction in the inlet will be high. 


Step 3. Measure the surface area of the bay at the mean water level, 
Ag, from charts or aerial photos. For most bays the surface area changes 
as the bay water level rises and falls because sections are flooded at 
high water levels. If the bay area change is significant, a bay area 
variation parameter, 8, is used to account for area of the bay, Abay » 
at any water level in the bay, hp, using the relation: 


Apay = Ao(1 + Bhp) , (2) 


where Ag is the bay surface area at datum, usually mean low water (MLW), 
mean sea level (MSL), or mean water level (MWL). 


Step 4. Specify the seawater level fluctuation as a function of time 
for the period of interest. Tide tables will give an estimate of the 
astronomical tide. Water levels can also be measured by a tide gage and 
stilling well (Seelig, 1977)%. Corps of Engineers and National Oceanic 
and Atmospheric Administration (NOAA) gages located at numerous points 
along the coast may also provide the desired water level information. 

In this computer program either the tide may be expressed as a sinusoidal 
wave with a period and amplitude or the levels may be described by instan- 
taneous sea level measurements at a constant sampling rate. 


Step 5. Determine the time step of input to the model for use in 
computations. As a lower limit, the time step, At, should be: 


_ =n 
At = Ved naa > (3) 


AX 


2SEELIG, W.N., "“Stilling Well Design for Accurate Water Level Measure- 
ment,'' TP 77-2, U.S. Army, Corps of Engineers, Coastal Engineering 
Research Center, Fort Belvoir, Va., Jan. 1977. 
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where L;, is the length of the inlet and d,,. is the maximum water 
depth in the inlet. A longer time step can be used for most tidal in- 
lets, and as an upper limit, the time step should be one-hundredth of 
the forcing wave period. 


Step 6. Document all input data using the computer format shown in 
the appendix. As a first estimate, set the flood and ebb entrance and 
exit loss coefficients to equal one (CDF = 1.0 and CDE = 1.0). Asa 
first approximation, Manning's n_ can be evaluated by the relation: 


me Cl S C2 ).. (4) 


where D is the local inlet stillwater depth. For depths greater than 

4 feet and less than 30 feet, Cl = 0.03777 and C2 = 0.000667; for depths 
less than 4 feet, C1 = 0.0550 and C2 = 0.005. The n for each grid may 
be different if C2 # 0. 


Step 7. For use with periodic forcing, run the program for several 
Sinusoidal cycles having the period and amplitude of the long wave of 
interest to approximate the hydraulic characteristics of the inlet-bay 
system. A Sinusoidal tide is specified in the model by giving the forc- 
ing period, T, in hours and the wave amplitude, A,, in feet, on card 
type 3 and by setting NPTS = 0 on card type 8 of input to the program 
INLET. Set ITABLE = 1 to obtain tables of instantaneous hydraulics at 
points throughout the water level cycle and set IPLOT = 1 to obtain a 
plot of predicted inlet velocities and discharge at sequential bay levels. 
These outputs will indicate the importance of the terms in the equation 
of motion describing water motion in the inlet. If temporal acceleration 
is small during most of the water level cycle, then startup transients 
will be small and the first or second cycle will contain little transient 
effect (NCYCLES = 1 or 2 in input data). However, if temporal accelera- 
tion is significant during more than 25 percent of the cycle, approxi- 
mately four cycles of model operation are required to eliminate startup 
transient effects (NCYCLES = 4). For aperiodic use such as with storm 
surges or rapidly varying wave size (e.g., tsunamis), run the model for 
the water level for approximately 10 hours before the time of interest 
to build up initial conditions in the model similar to the prototype. 


Step 8. Calibrate the computer model by varying Manning's n or 
flood- and ebb-loss coefficients. The seawater level fluctuation can be 
specified as a sinusoidal wave or in terms of an equal time series. For 
an equal time series, start and stop the series when the seawater level 
is at zero so that one or more complete cycles are described. Use at 
least 20 points to describe each cycle. The sampling interval in minutes, 
TDEL, and the number of points, NPTS, must be specified on card type 8 
and the water level data on card type 9. 


The model is calibrated using short periods of field observations by 
first comparing observed and predicted mean water velocities, if avail- 
able, at the minimum cross-sectional area region of the inlet. If the 
predicted velocities are higher or lower than observed, then the value 
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of n can be increased or decreased accordingly. When the computer 
model has been satisfactorily calibrated to predict inlet velocities, 
predicted bay water levels should be checked against measurements to 
assure that levels are being modeled correctly. If inlet velocities 
are not available, bay levels can be used to calibrate the model. 


Step 9. If additional prototype data are available, these data 
should be used to verify that the model adequately predicts inlet and 
bay hydraulics. 


Step 10. At this point the computer program is ready to use for 
prediction. Examples of the use of the computer program are presented 
in the following section. Input and output data, and computations are 
in U.S. Customary units. 


IIIT. EXAMPLES OF COMPUTER PROGRAM PREDICTION 


1. Cabin Point Creek, Virginia. 


Cabin Point Creek is a shallow natural tidal inlet that connects a 
bay to the lower Potomac River (Fig. 2) where the mean tidal range is 
approximately 1.5 feet. 


In this example, the model was calibrated with prototype river and 
bay levels and the calibrated model was then used to predict inlet veloc- 
ities, discharge, and bay level for a second inlet added to the system. 
The procedures for using the model are: 


(a) The inlet cross section was measured (Fig. 3) on 24 
May 1976, and is assumed to be representative of the 1,900- 
foot-long inlet. 


(b) The inlet is modeled using a grid system of three 
channels and two identical cross sections (Fig. 3) at either 
end of the inlet. 


(c) The bay area, A>» measured from a 74-minute U.S. 
Geological Survey (USGS) topographic map, was 3.5 x 10° 
square feet. For an increase in bay water elevation of 
0.25 foot, the bay surface area increases approximately 5 
percent because of marsh flooding. The bay area variation 
parameter, 8, can be determined from this information 
using equation (2), rearranged as: 


she. Abay i 
ae 3 (5) 


or, in this case, 


1 
= —— (1.05 - = 0.2 
p 0.25 ( 1) g 


Depth (ft) 


Bay Water Level 
Measurements 


Faoumen 2 


Figure 3. 
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Cabin Point Creek cross section. 
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(d) River water levels were measured at 30-minute intervals 
using a stilling well located near the inlet mouth (Fig. 2). 


(e) The time step was estimated as: 


1900 
At = === = 250 d 
30D OD seconds 


(f) Loss coefficients were specified as CDF = CDE = 1.0, and 
Manning's n was estimated as n = 0.055 - 0.005 D (recommended for 
depths less than 4 feet). 


(g) A preliminary computer run using a sinusoidal river tide 
showed that the inlet is controlled by friction effects and that 
temporal acceleration is not important. 


(h) The model was then run using the measured river water 
levels to force the model (Fig. 4). It was determined that the 
model adequately predicted bay levels. 


(1) No additional prototype data are available for verifi- 
cation of the model. 


(j) The model is now available to use for predictions of in- 
let hydraulics. In this example, a second inlet (inlet 2), is 
being considered for this site, so the model is used to predict 
hydraulics for the system with two inlets (Fig. 2). Procedures 
(a) and (b) are repeated for the second inlet. In this case, the 
second inlet is modeled by one channel and two cross sections so 
that the inlet has a length of 300 feet, a width of 50 feet, and 
a depth of 4 feet. These inlet data are put into the computer 
format, added to the program deck for the natural inlet, and re- 
run to predict conditions for the proposed two-inlet system. The 
numerical model predicts that addition of the second inlet would 
increase the tidal range and the tidal prism in the bay and would 
cause water velocities in inlet 1 to decrease (see Table). 


Table. Predicted Cabin Point Creek hydraulics. 


24 and 25 | Model prediction 
Tide May 1976 | for second inlet 


imnlec. i Inlet 1{ Inlet 2! 


Bay 0. 36 1.49 1.49 
(range in ft) 
Ebb -0.6 055) aio ¥ 


(maximum velocity in ft/s) 


Flood 0.9 0.3 ays 
maximum velocity in ft/s) 


1], = 300 feet, B = 50 feet, D = 4 feet. 


NOTE: Tidal range in the sea is 1.49 feet. 
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Cabin Point Creek sea and bay levels. 


2. Pentwater Inlet, Michigan. 


Pentwater Inlet is an example of a Great Lakes inlet controlled by 
vertical-walled jetties along the entire 2,000-foot channel (Fig. 5). 
Meteorologically generated seiches of Lake Michigan are the primary water 
level fluctuations causing reversing currents in the inlet. A model of 
Pentwater will be calibrated and used to estimate hydraulic response of 
the inlet to simultaneous lake seiching and river inflow. The procedures 
used in this modeling are: 


(a) A hydrographic survey of the inlet is used to describe 
the inlet geometry. 


(b) The inlet is modeled using one channel and six cross 
sections. 


(c) The bay surface area, measured from a hydrographic 
chart, is 1.81 x 107 square feet. The bay area does not 
change with bay water level because the bay has steep-sided 
slopes, so 8 = 0. 


(d) Lake Michigan water level measurements used to force 
the model were taken at 5-minute intervals on a tower located 
adjacent to Pentwater Inlet. 


(e) The model time step used is: 


2000 
A OCS ae 
t 359 x 15 90 seconds 


(f) Loss coefficients were specified as CDE = CDF = 1.0, 
and Manning's n was estimated by n = 0.03777 - 0.000667 D 
(recommended for depths greater than 4 feet and less than 30 
feet). 


(g) A preliminary run showed that temporal acceleration is 
an important term in the inlet equation of motion for Pentwater 
Inlet (Fig. 6). Therefore, several forcing cycles of model 
operation before the time of interest are necessary to eliminate 
transient terms due to startup conditions. 


(h) The model is calibrated by using Lake Michigan levels 
to force the model. An initial run showed that predicted bay 
level fluctuations adequately modeled observed levels (Fig. 7). 


(i) The model was not verified. 


(j) The model was used to predict inlet velocities, dis- 
charge, and bay levels for a 2-hour forcing wave with an 


IS 
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Approx. Scale (m) 
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Figure 5. Pentwater Inlet, Michigan. 
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Figure 6. Pentwater Inlet model prediction of monochromatic forcing 
(for a 2-hour wave with a 9.1-foot amplitude). 
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Figure 7. Pentwater Inlet model calibration. 


amplitude of 0.10 foot and a discharge into Pentwater Lake of 
2,800 cubic feet per second from the Pentwater River. The model 
predicted an average bay water surface elevation of 0.13 foot 
higher than the mean lake level, a bay water level fluctuation 
range of 0.25 foot, and a prism of water of 4.6 x 10© cubic feet 
caused by the seiche (Fig. 8). The inlet would always be in ebb 
flow due to river influence with a maximum velocity of -2.7 feet 
per second and a minimum velocity of -0.1 foot per second. Head, 
friction, and temporal and convective acceleration are important 
in the inlet equation of motion. 


IV. SUMMARY 


A computer program (INLET) based on_a numerical model (Seelig, Harris, 
and Herchenroder, in preparation, 1977)1 is presented for prediction of 
hydraulics where one or more inlets connect a bay to a sea. Two examples 
are given: (a) A tidal inlet forced by an astronomical tide where inlet 
channel friction is the dominant term in the equation of motion; and (b) 

a Great Lakes inlet with river inflow forced by lake seiching where head, 
friction, and temporal and convective accelerations are important at 
different points in the water level fluctuation cycle. The model can 
also be used for forcing other water level fluctuations, such as from 
storm surges or tsunamis. 


Another computer program (INLET2) is available for more complex 
systems of interconnected inlets, bays, and seas. INLET2 is an expanded 
version of INLET. Documentation and computer card decks for INLET2 are 
available from the Automatic Data Processing Division (CERDP), Coastal 
Engineering Research Center (CERC) . 


Details on model development and application, including additional 
examples, are reported by Seelig, Harris, and Herchenroder (in prepara- 
tion, 1977). 


1SEELIG, HARRIS, and HERCHENRODER, op. cit., p.- 7. 
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APPENDIX 
COMPUTER PROGRAM DOCUMENTATION (INLET) 
1. Program Description. 


The numerical model to predict inlet hydraulics is programed in 
FORTRAN for a CDC 6600 computer. The simultaneous differential equations 
are solved by a variable time step Runge-Kutta-Gill marching procedure. 
The organization of the computer program is shown in Figure A-1. A brief 
description of each routine follows: 


INLET is the main routine which controls input-output and calls 
subroutines to execute a specific task. Figure A-1 summarizes control 
throughout the program. The program is organized to accept up to three 
inlets connecting the bay to the sea, up to seven channels for each in- 
let, and up to eight cross sections (seven grids long). 


Subroutine HELM uses an iterative method of estimating the natural 
pumping period or Helmholtz period, Ty', for the inlet-bay system by 
neglecting friction in the inlet to give: 


Cig LY) A 
We S Zar Seon y Seats ebay 
BAe 
where L' is added inlet length due to radiation, and where L' is 
given by: 


Subroutine RKGS is a routine to solve simultaneous differential 
equations. This subroutine was adapted from the scientific subroutine 
package. 


Subroutine SETEQ evaluates the right-hand side of the equation of 
motion, one for each inlet, and the continuity equation between the inlet 
and bay for each step. This routine also evaluates the relative rank of 
the four terms in the equation of motion for flow in each inlet. 


Subroutine LEVEL determines the water level in the grids at each 
time step. The routine interpolates the level between the sea and bay 


based on the relative amount of friction in each grid cell. 


Subroutine TPWRTE writes hydraulic results from each time step on a 
tape or disc, so that this information can be used later by the output 
routines. 


Subroutine TABLE outputs a table of instantaneous hydraulics each 
time the routine is called. 
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Subroutine SEA determines the water level in the sea as a function 
of time either for a given sine wave or by interpolating equal time- 
series data. 


Subroutine WT1 determines the grid-weighting function by assuming 
that the flow is distributed across each section so that friction is 
minimized. This routine allows flow to cross channel boundaries, but 
assumes that this flow will be small, so the flow is neglected in the 
equation of motion. This weighting function is recommended for general 
use. 


Subroutine WT2 is similar to WT1, except that flow is not allowed to 
cross channel boundaries and that flow is distributed in each channel so 
that friction is minimized. 


Subroutine WT3 determines the weighting function so that flow is 
distributed equally in all grids. This is generally unrealistic, since 
it will be difficult to visually draw this grid system. However, this 
routine is useful since it provides an upper limit on frictional effects 
and therefore gives a lower limit of bay levels and inlet velocities. 
This weighting can be used to model simple geometry inlets where only 
one channel is used to represent the inlet. 


Subroutine CRIT prints a table of critical instantaneous hydraulics 
(i.e., at high water, low water, maximum velocity, and maximum discharge). 
This table is determined by storing a summary of conditions for each time 
step, then scanning this list for critical values. 


Subroutine GRPHC plots mean inlet hydraulics by scaling hydraulics 
in storage and plotting the time interval requested on a digital x-y pen 
plotter. 


Subroutine READIN is used by GRPHC to read data in storage and scale 
values for plotting. 


2. Program Input. 


The computer program (INLET) requires the following input of one deck 
for each inlet-bay system: 


Card Variables Format Description 
type 
1 ALABL1 4A10 first laa oF EMCO 
ALABL2 4A10 second line of title 
2 SLIO, Ze7O.55 Lilo 
NINLET number of inlets 
NCYCLES number of cycles 
IPLOT IPLOT = 1 for plot of results 


Ze) 


Card Variables Format Description 
type 
IWT weighting type 
IWT = 1 flow distributed to mini- 
mize (1 in card col. 40) 
ITABLE ITABLE = 1 for tables of instanta- 
neous hydraulics 
Gil, C2 Manning's n_ evaluated by: 
n = Gi = (€2 * Ds whexe Das steele 
water depth. If blank default 
values of C1 = 0.03777 and 
C2 = 0.000667 are assumed. 
ICONV ICONV = 1 (1 in card col. 80) 
3 3F10.5, E10.4, 
3F10.5, 2F5.1 
T forcing period (hours) 
DELT approximate time increment 
AO forcing wave amplitude (feet) 
AB bay area at datum (square feet) 
BETA bay area variation parameter 
ZETA inlet side slope D(z)/D(y) 
QINFLO bay inflow from sources other than 
the inlet (cubic feet per second) 
CDF an empirical flood-loss coefficient 
CDE an empirical ebb-loss coefficient 
4 Ziel OR EALORO 
IC number of channels 
IS number of cross sections 
QINT estimated inlet discharge at the 
time the model starts 
5 (one card 10X, 7F10.5 
per 
section) 
At cell cross-sectional areas at the 
ends of each cell at datum (square 
IEE) (SEQ Wales (N52) 
6 (one card 10X, 7F10.5 
per 
section) 
B! grid cell widths for the end of 


each cell (feet) (see Fig. A-2) 
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Card Variables 


Format Description 


7 (one less 
card than 
sections) 


IL Y 


IOXS TIO >S 


lengths of the sides of cells (see 
Fig. A-2) (one less card than num- 
ber of sections; one more value per 
card than the number of channels) 


For card types 5 to 7, there will be one card for each cross section 


of the inlet. 


The first card will be for the first cross section, i.e., 


the section closest to the sea, and the last section is adjacent to the 
bay. The first value on each card will correspond to the first channel 
adjacent to land; the last value on each card will correspond te the 
last channel also adjacent to land (Fig. A-2). 


For more than one inlet connecting the bay to the sea, repeat card 
types 3 to 7 for each additional inlet. 


Card Variables Format Description 
type = 
SD NeilL 34X, F6.2 water level sampling interval 
(minute) 
NPTS OX, US number of sample points = 0 for no 
data 


9 (optional--no cards if NPTS = 0 from card type 8) 


YG 


eight water level values per card, 
as many cards to include NPTS points; 
start the model at a time when the 
sea level is zero. Use 25 or more 
points per forcing cycle for best 
results; i.e., levels at 30- or 15- 
minute intervals for a 12-hour tide. 


10 (optional--two plot cards, first card used only if 
IPLOT = 1 on card type 1) 


XO 
XF 


SCALX 
YLO 
YL 


YLSCAL 


SITs 5 5/ sH21055, 0) 
starting time of plot (hours) 
ending time of plot (hours) 
time scale (hours per inch) 
minimum value of water levels (feet) 
overall height of plot (inches) 


scale of water level height (feet 
per inch) 
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Card Variables Format Description 
type 


YRO minimum flows (thousand cubic feet 
per second) 
YRSCAL scale of flows (thousand cubic feet 


per second per inch) 


Second card 


YVO minimum velocity (feet per second) 

YVSCAL scale of velocities (feet per 
second per inch) 

SCALE scale factor for total plot size 

IQ IQ = 0 for no plot of inlet 
discharge 


11 If a plot is requested, repeat card types 8 and 9 for observed bay 
levels to compare with predictions (card type 8 required; use 
NPTS = 0 for no observed bay levels). Only one set of card types 
10 and 11 will be required for plotting even though the system 
modeled may have more than one inlet. 


12 End of file card. 


The inlet data for a computer run of Masonboro Inlet, North Carolina, 
are shown in Figure A-3. 


3. Program Output. 


The types of output include: (a) A summary table of grid dimensions, 
input parameters, and the Helmholtz period of the system estimated assum- 
ing there is no friction in the inlet; (b) (optional) summary tables of 
instantaneous inlet hydraulics; (c) (optional) a pen plot of inlet hydrau- 
lics; and (d) a table summarizing critical points throughout model opera- 
tion, such as high water, low water, point of maximum discharge, and 
maximum velocity. Samples of input and output for the Masonboro Inlet 
run are given in Figures A-4, A-5, and A-6. 


4. Computer Program. 


A listing of the computer program (INLET) follows the sample output. 
The program was written in FORTRAN IV for a CDC 6600 computer with plotter. 
Control cards, plotting instructions, and file controls may have to be 
changed for other computers. If no plotter is available, the subroutine 
GRPHC and the call to the subroutine in the main program may be removed. 
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MASONKBORO 10469 
CoFa2, 
1 2 
2500 200. 2.35 e20000F4N9 
7 ©20000, 
AY PERRI 55106 U5706 
42 9725, THRS, SABO, 
Ag BBA. SAA, SoPs.e 
Wg 940, 2525. 19030. 
AS EN. VORA, 5079, 
46 3770, SASSO, 5540, 
47 Ute O41%, KUN» 
By BAe AO. 200, 
82 132%. 14006 310.6 
BS ENO. 1SRUe 20. 
Bu 359. 430. 450. 
85 PRO 1S. Aen, 
Be RUD. AON, UPne 
87 OA0. 6706 470. 
ui ASN. 900, yoann, 
L2 75%6 950. 109, 
L3 nHOe 559. Qn. 
La 510. 700. Asn, 
LS une BOO. 950. 
Lo Dann, 2100. 2'04, 
GAGEN 9712/49 “ASONBIRN OFLTs 30, 
21.39 o140 e1,h45 =1.60 
0034 ne 1.?9 1.70 
264) 2.6ce 1.9} Gia) 
=0,98 e132 =1.55 =1.62 
20020 e36 0293 1.40 
204A 2.09 1.97 1056 
=0.9 °1,3 
0. 22. 26 «3. 
she eo le 
NO BAY 
EoR 


Figure A-3. 


euro, 
e1une 
3790, 
SPAS, 
“nen, 
395, 
GNDMs 
90. 
190. 
2400 
sun, 
35%. 
U0. 
e606 
1900, 
on, 
1950, 
gon. 
60%. 
3604, 
NUM 50 
o1l.3R 
200A 
lo 
"1.60 
1o74 
lel 


6. 


1900, 
N00, 
1300, 
900. 
20%, 
3400. 


=0,98 


2.53 
0,50 
oldu 


2010 
006 


Jo 


Sample of input data for a computer run of 
Masonboro Inlet, North Carolina. 
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MASONBORO 1969 


YES? 


CONTROL CARDS 


1 1 0 2 
25.00000 299.00000 201500 ce0vne+o9 


SUMMARY OF INLET GRIND CHARACTESISTICS 


a 
SECTION 1 
CHANNEL & { 2 
ARFACFT>) 1900205 6697.5 
wIOTHC(F Tr) 2160.0 1040.0 
DEPTH(Fr) 8,Aa0 6,44 
LENCFT) 875.0 95060 
N 003149 00335 
SECTION 2 
CHANNEL 8 { ? 
AREA(F TD) 6402.5 6767.5 
WIOTH(F Tr). 910.00 13900n 
DEPTH(FT) 7,04 4,47 
LENCFT) 850.0 975,90 
N ©0331 003458 
SECTION 3 
CHANNEL 3 1 2 
AREACFT2) 201000 4967.5 
WINTH(FT) 42500 2905.0 
DeEPTN(FT) 4,73 4,82 
LEN(FT) 495.0 72500 
N ©0346 00348 
SECTION u 
CHANNEL 3 1 2 
AREAC FT|) 72000 2780.5 
WIDTM(FT) 31500 29960 
DEPTH(Fry 2029 9.59 
LEN(FT) 600.0 77509 
N 00362 04314 
SECTION 5 
CHANNEL @ { ? 
AREACFT >) 213500 4443.0 
WIOTH(F Tr) 560.0 520.00 
DEPTH(FT) 3,81 4,54 
LENCFT) 600.4 AIS.N 
N 00352 ©0321 
SECTION 6 
CHANNEL o@ { ? 
ARFACFT>) 408000 OP 3000 
wWIDTHCFT) 910.0 780.0 
DEPTHCFET) 4,48 7.99 
LENCFT) 2350.9 2100.9 
N 00348 00324 
FORCING PERyONS 25,00 HOURS 
THELM(APPROy) a 3,17 HOURS 
TF/THe 7,84 
INLET LENGTH ADNED LENGTH 
1 goe20e5 170964 
TOELo MyNe 30.00 NPTSz 54 
21,39 @1,60 21.65 91,40 -1,58 
Pod) 2,22 1691 1.650 1,00 
2,20 236 oO Tot Vor as 
2.970 ef 40 
Figure A-4. 


INLET NUMBER 4 
6 


1 0.00000 0.00000 
220000 001330 0.00000 220 0,0 
3 4 
5125.0 228000 
285,0 9500 
17098 24.00 
1000.9 100000 
2025A 00218 
3 a 
5452.5 2920060 
295.0 180.00 
19,16 16,232 
1900.0 10000 
2025A 00269 
3 a 
7827.5 449265 
365.0 40000 
21,45 11.23 
975.0 112509 
20235 00303 
3 a 
7554.5 468265 
365.0 44u5en 
20070 10.52 
A75.0 900.0 
20240 00308 
5 u 
$204.5 400205 
35n.0 40500 
14,47 9,A8 
77500 4004 
20279 ©0312 
3 4 
665.0 3962.5 
545.0 30004 
12060 11.401 
2850.9 3500.0 
2029u s0304G 
98 6.49 = 98 osu eRe 1.29 $,.70 2.98 2,33 2.48 2,50 
SQ MeN we eSU 098 HO 1432 M1055 ef ebe @1,60 a1 uu 1403 @,69 
0 Pod BohY Pode 2529 VoD? Dod Vole 060 010 #40 


Sample output from INLET (summary table for 
Masonboro Inlet input data). 


2g 


TES 


CRANKEL 


CO ie ie ee a Dl el ll cee 


4 
TREMP 


MEAN VELOCITY AT THE 


e 4AuURS 


ISLET 


= 6.000 DELTs SEC s 


1 


SEa LEVELeFT= 2.08 
Sav LEVELeFT= 123 
DISCHARGE*CFS= ,SuslETOS 


206 


2.08 
033 
2802, 
005 
2000 


2002 
093 
8993, 
16 
001 


2090 
4.94 
3125R, 
057 
093 


200d 
3050 
11772, 
021 
04 


Bay AREAS ,2UG3E+09 FT2 
SECTION 1 
FRI¢ 004 
LEVEL 2.08 
VCFPS) e012 
G(CFS) 2B02. 
wETGRT 005 
FRI¢ 290 
LEVEL 2.96 
vc(Fes) 1201 
Q(CFS) 899%, 
WEIGHT 016 
Fre 001 
LEVEL 2.096 
v(FoS) 5.40 
Q(CFS) 31238. 
WEIGHT 057 
FRIC 0035 
LEVEL 2.07 
Vc(FpS) 4260 
9(CFS) 117726 
WEIGHT 02k 
FRIC A) 
4Ccs e6 CONV ACCS 


400.00 


3 
007 


2.46 
204 
2802. 
20S 
000 


1,94 
1.52 
&99%, 
o!6 
202 


1,95 
3.63 
T125R, 
057 
002 


fe 


a 
042 


1070 
2014 
2R02. 
09S 
010 


{066 
Poll 
8993, 
016 
e010 


1083 
Wo77 
31236, 
057 
ell 


1075 
2013 
11772. 
oe} 
010 


32,4 HEADS 219000 FRICS 


MINIMUM AREA SECTIONS 


2.97 FT/SEC 


eel 
091 
67.0 


A4TNS 


6 
031 


1.26 
053 
2802. 
00S 
201 


1,29 
1,2u 
8993, 
216 
203 


1,42 
4,07 
312386, 
057 
220 


1,37 
2.62 
11772. 
eek 
208 


18.29.73 Fre 


FRICTION 


el2 


023 


Figure A-5. Sample output from INLET (summary table of instantaneous 


hydraulics for Masonboro after 6 hours of model time). 
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Figure A-6. 


245111 
24.223 
25.2000 


SUMMARY TABLE OF MYDRAULICS 
INFLOW 


HS 
FT 


KCFS 


0.000 
02000 
02090 
02000 
02000 
02000 
02000 
0-000 
02000 
0000 
06000 
00000 
02000 
00060 
00000 
00000 
02000 
020900 
02900 
02090 
00000 
00090 
0°V00 
02090 
00000 
0°00 
0000 
00000 
02000 
02000 
40000 
02000 
02000 
02000 
0090 
02000 
90000 
02000 
02009 
0°000 
92000 
00090 
02000 
G2000 
02090 
02000 
0000 
02000 
02090 
02000 
00000 
02000 
02000 
02000 
92000 
02000 
0000 
00000 
00000 
02000 
02000 
00000 
02000 
02000 
02000 
02000 


* CRITICAL POINT VALUE 


Sample output from INLET (table of critical points for 


the model time: high water, low water, etc., for 


Masonboro Inlet). 


3| 


HA 
rT 


=.239 
=,951 
1.562% 
Ong) 
2,456 
2516 
2698 
.788 
2876 
2967 
1,056 
1,145 
1,234 
treesies 
e.147 
2,46e" 
1.191 


VEL 
FPS 


=3,86j1% 
=20919 
0053 
2ce463% 
2eN81* 
20920” 
20940" 
26945m 
Oe 94he 
20957 
209hRe 
2e976e 
2o.9?ue 
2c95Au 
2e1Su 
2086 
=3.308 
25.53 7e 
23.3628 
o§,3ADu 
=3.39A" 
=3.U)ju 
23,4eon 
23,429" 
=3,4s ane 
23,43 sn 
23,430" 
5.427% 
230420" 
=o3.403" 
210764 
220738 
10880" 
2e9IUm 
30020% 
3o036" 
3B.0RGe 
3.002% 
32004 
Be053e 
22994 
20973 
20965% 
20967« 
20969" 
20942 
20267 
2.016 
=20.904% 
2209? {x 
220942 
=3e3IUH 
230415" 
=3.429% 
23.440" 
=5o.44Gn 
=30456mu 
=3o459" 
=3.46j«% 
=3.46{* 
23.463" 
=504b6o" 
@3,.45um" 
=30435% 
3,409" 
22.599 


INLET 


@55,1460* 
239,568 
2685 
347,947 
%B.6351 
50.c6 
51,646 
52.193 
52.090 
53.252 
53.684 
54,u4} 
5u,806 
64,6098 
u1.977 
1.734 
@55,734% 
°55.713 
255.6007 
255.405 
-S5.177 
254,870 
254,519 
e54,12o 
=53,0580 
53.170 
=52,600 
252.037 
51,412 
250.697 
222.798 
-.923 
25.949 
50.979 
52.006 
52,865 
53,6808 
53.085 
53,720* 
53.719 
5y,44e2™ 
53.468 
53,749 
54.204 
54,048 
54, 883% 
44.165 
=,3ie 
252,026*% 
252.545 
=52.477% 
256,659% 
256,474 
256.194 
255.8356 
255,460 
255,044 
e54,588 
254.092 
253,574 
253.063 
252.5156 
251.870 
251.065 
=50.167 
235.948 


Listing of the computer program INLET. 


PROGRAM INLETCINPUT OUTPUT es TAPESSINPUTyTAPEG=OUTPUTeTAPEQeTAPEXO® INLET 


1 TAPE3sPUNCHETAPES) INLET 

C PROGRAM NUMRER 720K6R16S0 (INLET) ANALYSES AND PREDICTS INSTANEFOUS INL INLET 
C RYDRAVLICS USING A LUMPED PARAMETER SCHEME (SEE SFELIGs HARRIS AND INLET 
C HERCHENRODERe 19765 [A GENERALIZED LUMPED PARAMETER MODEL OF INLET INLET 
C HYDRAULICS[+ A ORAFT CERC REPORT) INLET 
REAL LeLENGTHeLINeLXoNoNX INLET 
COMMON/NUMS/NT eo GoNTINLETOICH(3) ep ISE(3) eQROL(7O7)oBKTOTIODE 797) 0 INLET 

1 AC7T¢7)eNO707) 00797) OVE 707) 9Q0797) oHSeHBoH( 797) 0 1CoTSoAMINI (3) 0 INLET 
1BMINT C3) eLINe 9X03) eAINFLUCARAYOLENGTH( 3) INLET 
COMMON/NUML/Y(5) eDERY(S) oXONTOIWT eo ZETAGHH INLET 
COMMON/NUM2/BX (30797) oDXC3Se70 7) OHX(SeTo PF) OWNC Se 70O7IOLX(S07e7)ONK(S INLET 
19797) INLET 
COMMON /NUMS/A0eTOARG RETA INLET 
COMMON/NUMU/RNK (394) INLET 
OIKFNSION CORL(3) INLET 
DIMENSTON ALABL1(4) eALABL2(4) o IBUF (1000) oNUMBER(20) INLET 

3370 CONTINUE INLET 
bO 2193 Ilete3 INLET 

2195 QX(TI1) 31. INLET 
c Ge ACCELERATION OF GRAVITY INLET 
Gs32,2 INLET 

OO 1241 151020 INLET 

1211 NUMNAER(T)eY INLET 
WRITE (602937) INLET 

2937 FORMATO //04 X90 (pocteareceseccessccceces () INLET 
REAN(Se1167) C(CALARLI(T) 018194) TNLET 
REAN(5e1167) (CALABLA(1)°181,4) INLET 

1167 FORMAT(UAIN) INLET 
WRITE (601168) (CALARLICI) 912194) INLET 

WRITE (60116A) CALABI2(1) 918104) INLET 

1168 FORMAT(AXe4A10) INLET 
WRITE (601268) INLET 

1268 FORMAT(seSxe (CONTROL CARDS () INLET 
C READ CONTROL CARDS INLET 
c INLET 
REAN(S94014) NINLETeNCYCLESs IPLOTeIWToITABLE OC Ce INLET 

WRITE (601942) NINLETyNCYCLESy IPLOTy IWTeITABLEeCleCe2 INLET 

1011 FORMAT(ST1% 92F10,5) INLET 
1012 FORMAT(4{xXeSI10,2F 10.5) INLET 
C NINLET2THE NUMBER NF INLETS INLET 
C NCYCLES= NUMRER OF TIDAL CYCLES INLET 
C IPLOT (1 FGR A PLOT OF MEAN HYDRAULICSe 0 FOR NO PLOT) INLET 
C IWT IS A PARAMETER NESCRIABING THE TYPE OF WEIGHTING DESIRED INLET 
C IWTel! FOR FLOW WEIGHTING TO ACHIEVE MINIMUM FRICTION INLET 
C IwTs2 FOR WEIGHTING FoR MINIMUM FRICTION WITH NO FLOW ACROSS CHANNELS INLET 
C IWT=3 FOR FOUAL FLOW TN aLL GRIDS TO GIVE MAXIMUM FRICTION INLET 
C ITABLE&;) FOR A TABLF OF OLITPUT INLET 
C CleC2 NECieCe * D, IF C1 AND C2 ARE ZERO THE MASCH VALUES OF INLET 
C Ci £.03777 AND C2e.,090K67 ARE USED INLET 
TF (C1 0E9e0,0eANDoC2.EN000) C2B 0.000667 INLET 
IF(C10EQ00,) C180003777 INLET 

Cc INLET 


S12 


ee ee 
WEwrw COO ONoNMEWMN 


16 


1 FORMAT(AII0) 
REAN(Set11) TeDFELT, AO gAHOHETAGZE TAS QINFLO 
VRITECOolLIE) TeMELTe AN ABO KETAGZETAGGINFLO 
LL}  FORMATCBF LO eS0E10ed dF 1005) 
TSTIDAL PFEKTODe HRS (LATER CONVERTED TO SECONDS) 
DELTZSESTIMATED TIME STEP SEC 
Ans SFA TTOAL AMPLITUNE GFT 
ARe BDAY AREA AT THE DATIIMe SQUARE FEET 
BETAS AAY AREA VARIATTON PARAMETER C€ N(AB)/0(HB)) 
ZETAS CHANNEL SLOPE (N(¥)/0(K)) 
WINFLOS INFLOW INTO THE BAY FRUM OTHER SOURCES (FT3/SEC) 


ananannnn 


ENDeT*®NCYCLES* 3600, 
IF(7ETA,LE,O,)ZETABI E25 
NT80 


c 
C READ IN INFORMATION OF FACH INLET 
DO 4410 NIsteNINLET 
IUNTT=B4NI 
REWARD TUNTT 
REAN(5e1) ICoIS 
C Ce NuMBER OF CHANNELS 
C Se NUMBER OF INLET CRNSSeSECTIONS 
IFC TCoGTe7?,YReISeGT,7) WRITE (601671) 
1671 FOKNATC///¢5X0 (#e#en TOO MANY GRIDS FOR DIMENSIONS ty//) 
ICH(NI) IC 
C READ SFCTION AREAS ( ONE CARD PER SECTION) 
PO § Te1els 
5 REAN(Se2) (ACT) J) eJete IC) 
2 FORMAT(10X%97F10,5) 


C 
C READ SECTION WIDTHS CONE CAKD PER SECTION) 
00 & IaieIS 
6 REAN(5e2) (CHCTed) oJaioIC) 


C 
IcPysICe]} 
TSMislSe1 
C READ LENGTHS (ONE MORE LENGTH PER CARD THAN CHANNELS) 
C ( ONE LESS CaRD THAN THE NUMBER OF SECTIONS) 


00 7 IeyelISmji 
7 REAN( Se?) (LC Led) oJeteICPh) 


C 
C INITIALTZE VARYABLES TO BFGIN ITERATION 
C NUMHER MF GRINS ALONG THE CHANNEL IS ONE LESS THAN THE NUMBER OF 
C CRUSS@*SECTIONS 
66 1Sé1s91 
TSECNI)3IS 
ISHyelSe} 
WRITE (603678) NI 5 
3O7B8 FORKMAT( /95Xe (SUMMARY OF INLET GRID CHARACTERISTICS (o/ 
1 15x¥0 (INLET NUMBER (913) 
wWRITE(601) ICe1S 
DO 10 I#ieJs 


33 


INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
TNLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 


M1 
1297 
leeli 


HO 4) Jeresc 

LENGTH(NI) SLENGTHONT) +L 1ed)/FLOAT(IC) 
ACTeJIE(ACT oJ) +ACTe1eJ) 20 
LOLeJIECL CT eS) +l Closet) 7206 

HCle JISC HC Te SIHACTetod) I/e0 

OCT eJISACT eI) KCI 0) 

NCI eJ)eCi=Ce*0(Ted) 
LXCNTeTeJ) ol (Jed) 
RXONTe Te J) eA I 9 J) 
DXCNToTe J) eC Ted) 
NK(NToleJ BNC TeJ) 
4XCNToTeoJ) a1 e/FLOAT(IC) 

CONTINUF 

*RITE(CO01297) I 

FORMAT(/e1xXe (SECTION (013) 

wRITEC 691221) (NUNHFROIT) ol TetelIC) 
FORMAT(SXe (CHANNEL s[(el0I109/) 


C PRINT A SUMMARY TABLE OF GEOMETRIES 


1971 
1972 
1973 
1974 
1975 
10 


WRITE (601971) (ACTe J) eJELeI“) 
WRITE C6O01972) (ACT eI) oUBl1eIC) 
YWRITECGO61975) (NC Ie J) edFlelc) 
WRITE(601974) (LO19.1) oJELeIC) 
®RITEC601975) (NCIe Js) oJBle1C) 
FORMAT(SXe (AREACFT2) (o19F 10,1) 
FORMAT(SXe (WIDTHCFT) (010F 10,1) 
FORMAT(SxXe (DEPTH(FT) (o1Xel0F1002) 
FORMAT(SxXe (LEN(FT) (oeXol0FI0e1) 
FORMAT(SXe (N[o1OXe10F 1004) 
CONTINUE 


C FIND AREA AND *IOTH AT THE MINIMUM SECTION 


109 
1110 
c EST 


201 


1337 


2268 


2269 


AMINICNI) 299, Eet2 

00 109 Jei,I8 

AAgg, 

BRS0, 

00 108 Jet,IC 

AAMAAPA(TI oJ) 

BRERBSH (Tod) 

IFCAAeGT,AMINI(NI)) GO TO 109 

AMINICNT)BAA 

BMINICNI) SKB 

CONTINUE 

CONTINUE 

IMATE THE INLET*=RAY HELMHOLTZ PERIOD 

CALL HELM(THELMy AByCORL) 

THTFST/THELM 

WRITE (69201) Te THELMe THIF 

FORMAT(1Xe [FORCING PERIODE [gf 7.20 HUURS lo 
L/olXe (THELMCAPPROX)a(0FB,20[ HOURS (0/ 

1 1Xe (TF/THS Col 0K9F 6.2) 

WRITE(601337) CCJ*LENGTH(J) eCORL(J)) oJB1 oNINLET) 


FORMAT( «3X9 fINLET LENGTH ADDED LENGTH(s (/o4KolTeolXe 


1 FO,192XeF6e1)) 

TeT*3o00, 

CALL RKGSCEND*ODELToNINLETeQTNFLOoITABLE eT) 
DELTEEND/FLOAT(NT) 

DO 2269 NIsieNINLET 

HHEHS 

MRITEC6Oe226R) NI 

FORMAT(//010Xe (SUMMARY TABLE OF HYDRAULICS INLET (915) 
IUNTTENTI4R 

CALL CRITONTeDELTSTUNIT9TeNCYCLES) 
IFCTPLOTSEQs1eANDeNT,FQol) CALL PLOTS(IAUF 9100003) 
IFCTPLOT,EG,1) CALL GRPHCCALABL1eALABL2¢DELToTUNITONI) 
TFC TPLOTsEQe1eANDeNToFQeNINLET) CALL PLOT(0.9059999) 
CONTINUE 

STOP 

END 
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INLET 
INL&T 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 


108 
109 
110 
111 

112 
113 
Lia 
TS 
116 
117 
118 
119 
120 
lei 
lee 
123 
lea 
125 
126 
127 
128 
129 
150 
131 
132 
133 
134 
135 
156 
137 
158 
159 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
151 
152 
153 
154 
155 
156 
197 
158 
159 
160 


16} 
lo2 
1o3 
lou 
165 
166 
167 
168 
169 
170 
171 
172 
175 
174 


C ROUTINE TO SOLVE A SET MF SIMUTANEOUS DIFFERENTIAL FQUATIONS 
C ADAPTED FROM SCIENTIFIC SUBROUTINE PACKAGE s 


li2e 
1 
2 
3 
4 
5 
6 


SUBROUTTNE RKGSCENDeDELTeNINLETeQINFLO,TTABLE oT) 


IBMe 


COMMON/NUMI/¥ (5) oDERY(S) oXONTOIWToZETAGHS 


COMMON/NUMU/RNK (Sed) 


DIMENSION AUX (B95) 9A(8) 0808) 9C(8) oPRMT(5) pAMINIC3) 


NOIMENINLET#Y 

PRMT(1) @ Jeo 
PRMT(2)=END 
PRMT(3)EDELT 

PRMNT(4) & 43 

TFC T.GT, 360900) DELTHE 36006 
IFCT,LE,36000,) DFLTBBT/9%. 
NO 1122 JNzloNINLET 
YOUN) 70,01 
DERY¥(JN)E0,901 
Y(NOIM)e0, 
DERYCNDIM) a! O@FLOAT(NINLET)*0,001 
DO 4 ITzteNDIM 
AUX(B0T)80,0666666678DERY (1) 
kfPRMT(}) 
XENNEPRMT (2) 
HEPRMT(3) 

PRNT(5S)300 

CALL SFTEQ( AMIN) 

IF CHR CXEND@X))3A0 S702 
CONTINUE 

A(1)29.5 
A(2)=0.2928932 

AC S$)21.707)07 
4(4)=9,16666667 
B(1)se. 

Bl2yal. 

B(Sysi, 

b(4)s2, 

C(1)20.5 
C(2)=0,2928932 
CC3)sl, 707107 
C(4)30.5 

DO 3 ITsteNnDIM 
AUX(19T)SY(1) 
AUX(2eT)EDERY(T) 
AUX(3eT)B0, 
AUX(601)80, 

IREC30 

HEHeH 

IHL Fee} 

ISTFEPEO 

IF Nns0 

CONTINUF 

TFC CXtHeXEND) #H) 70605 
CONTINUE 

CONTINUE 
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1970 


INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
TNLET 
INLET 
INLET 
ENLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
TNLET 
INLET 
INLET 


175 
176 
177 
176 
179 
160 
161 
162 
163 
164 
165 
166 
167 
168 
189 
190 
191 
192 
193 
194 
SS 
196 
197 
196 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
eid 
eS 
216 
217 
218 
219 
220 
221 
222 
e223 
224 
225 
226 
2e7 


17 


19 


18 


20 


21 


HEXFNDeYX 
IENps! 

CONTINUE 

CALL SEA(HSeX) 

CALL TPWRTECNINLETe Xe HSOQINFLOD Yo AMINI eRNK ONT) 
IFLAGIExX/DELTB 
IFCIFLAGL oNEgIFLAGS,ANDeITABLEoEGes) CALL TABLE 
IFLaGesTFLaGi 

IF CPRMT(S))4008,40 
CONTINUE 

ITEST=0 

CONTINUE 
ISTFPEISTEPSS 

Jal 

CONTINUE 

AJSA(J) 

RBJ=A(J) 

cJer(J) 

0O 11 Is{enDIM 
RisHeDERY(T) 
R2SAT*(RleBIeAUN(607)) 
YC I)ev(t)¢re2 
R2ER2¢KD+RD 
AUX(Cor1) SAUX (6eT) #RAeCUFRI 
TF (J24)12015015 
CONTINUE 

JoJe{ 

IFC J=3) 13014913 
CONTINUE 

X2X4005*H 

CONTINUE 

C4LL SETEQ(AMINT) 
GO TO 19 

CONTINUE 

IFC TTEST) 16016920 
CONTINUE 

00 17 IzlyNOIM 
AUX(UeT)sY(I) 
ITESTE4 
ISTEPZISTEP+ISTFPa2 
CONTINUE 
IHKLFSIHLF +] 

XBXoh 

H20,5*H 

DO 19 IsieNNI™ 
YCI)zAuX(1,T) 
DERY(])szAUx(2e1) 
AUX( 691) ZAUX(39T) 
GO TO 9 

CONTINUF 
IMODSEISTEPs2 

IF (ISTEPeIMODSIMOD) 24,2302) 
CONTINUE 
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INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 


228 
2e9 
230 
25) 
232 
233 
234 
235 
256 
237 
238 
259 
240 
euy 
e4e 
243 
244 
245 
246 
2u7 
248 
249 
250 
251 
2o2 
253 
254 
255 
256 
257 
258 
259 
260 
26) 
262 
263 
eo 
265 
266 
267 
268 
269 
270 
e714 
272 
273 
274 
275 
276 
e77 
278 
279 
280 


22 


eu 


e7 


29 


31 


40 


23 


25 


26 


26 


30 


32 


33 


34 


35 


36 


37 


36 


39 


CALL SETEQ(AMINT) 

O09 P2 IaleNDIM 
AUX(SeI)8Y(1) 
AUX(7TOT)SDERYCT) 

GO TO 9 

CONTINUE 

OELTSO, 

DO Au TsleNDIM 

DELTEDEL THAUX (Bol) #ABS(AUX(UeT)@¥(1)) 
IF (MELTSPRMT(4)) 28928029 
CONTINUE 

IF CTHLF E10) 26036056 
CONTINUE 

ODO A7 TelenDIM 
AUX(4eT)3AUX(SeT) 
ISTFPZISTEP+ISTFPOY 
XEXeH 

ITENn&0 

GO TU 1A 

CONTINUE 

CALL SETENCANINT) 

00 29 IslenDIM 
AUXCLOT)SEY(T) 
AUX(2e1)SNERY(T) 
AUX(SeT)SAUX(OeT) 
Y(TysAux(Se1) 
OERY(I)sAux(7eT) 

CALI SEA(HSeX=H) 

CALL TPWRTECNINLETe XeHeHSeOTNFLOoVeAMINTORNKONT) 
IFLaG1S(XeH) JNELTB 

IFC TFLAGI,NEIFLAG2,ANDoITABLE,EG,1) CALL TABLE 
IFLAG2=IFLaAGl 

IF (PEMT(S))40030040 
CONTINUE 

OO 31 IslenhI™ 

YC1)sAux( tel) 
DERY(I)sAUxC20T) 
IRECSIHLE 

IFC TEND) 32932039 
CONTINUE 

IHLFeIHl Fey 
ISTFPSISTEP/2 

HZHeH 

IF CTRLF I 4033933 

CONTINUE 

IMUNZISTEPs2 
IFCTSTEP@IMOD@IMOD) Ug ued 
CONTINUE 

IF (NEL T©0,02"PRMT(4)) 3503504 
CONTINUE 

ITHLFSIHLF®} 
ISTFPZISTEP/2 

HEASH 

GO mT 4 

CONTINUE 

THLFe11 

CALL SETEQ(AMINY) 

GO TO 39 

CONTINUE 

THLF ele 

GO To $9 

CONTINUE 

THER EIS 

CONTINUE 

CALL SEACHSeX) 

CALL TPWRTIECNINLE Te X¥eHSOQINFLOeVoAMINI go RNK ONT) 
IFLAGIEX/OFLTB 

IFCTPLAGL NE IFLAGZ ,ANDeITABLE,EG,1) CALL TABLE 
IFLAGCBTFLAGI 

CONTINUE 

RETURN 

END 
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INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INL! T 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
TNLET 
INLET 
INLET 
INLET 
INLET 


281 

262 
283 
2bu 
285 
286 
267 
268 
289 
290 
291 
292 
293 
29u 
295 
296 
297 
298 
299 
300 
301 
$02 
303 
304 
305 
306 
307 
308 
309 
310 
311 
312 
313 
314 
315 
516 
317 
318 
319 
320 
321 
3e2 
323 
324 
3e5 
326 
327 
328 
329 
330 
331 
332 
333 
354 
335 
336 
337 
338 
339 
340 
Jay 

342 
343 
344 
$45 
346 
347 
$48 
349 
350 
351 

352 


SUBROUTTNE SETEQCAMIN)Y INLET 353 


C RNUTINE TO SEUP THE EQUATIONS FoR THE RIGHT HAND SIDE OF THE EQUATIONS INLET 354 
C MOTION AND TO DFTERMINE THE RANK OF THE TERMS IN THE EQUATION OF MOTIO INLET 355 
REAL LeLENGTHeLTNeLXoNoNXeLE INLET 356 
COMMON/NUMS/NT oe GeNINLETOICH( 3) 9 ISE( 3) eQReL (707) 0B( 707) 9D( 707) 0 INLET 357 

1 AC7e7TI ONC 707) 06797) 0 V6797) 900797) oHSeHBOH( Zo 7) oICoTSeAMINI(3) 0 INLET 358 
LBMINTOC3)eLINeQX 05) eQINFLOeARAY ep LENGTH(3) INLET 359 
COMNMUN/NUMI/Y0(5) eDERY(5S) oXONTO INT oO ZETAQHH INLET 360 
COMMON/NUM2/BX (39797) oOXC 30707) OHK( S07 07) OWK(So7O7) OLX(Se7To 7) ONX(S INLET 3o1 
19707) INLET 362 
COMMON /NUM3/A0,TOARY BETA INLET 363 
COMMON/NUMU/RNK (304) INLET 3ou 
DIMENSION AMIN(3) INLET 365 
G=32,e¢ INLET $606 

DO 220 NIB} 03 INLET 367 

NO 119 TH1—94 INLET 368 

119 RNKONToT) 30. INLET 369 
220 CONTINUE INLET 370 
CALL SEACHS eX) INLET 371 
HHEHS INLET 372 

C FIND THF BAY AREA INLET B7/s} 
HRSYv(NINLFT4+1) INLET 374 
ARAywABY* (1 ,¢BE TAHA) INLET 375 
QTso, INLET 376 

C SET UO FEUATIONS FOR EACH INLET INLET 377 
O09 100 NISyeNINLET INLET 378 
AMIN( NT) 2599999999999, INLET 379 
Gnsy(NT) INLET 3560 
QTENT+09 INLET 381 
ICEICH(NI) INLET 382 
IS=TSE(NI) INLET 563 
LEE, INLET 364 

00 95 I2leIS INLET 385 

DO 94 JaleIC INLET 386 

NOT e J) SNKONT oT 9 J) INLET 387 

LOT es JIELXONToT oJ) INLET 368 
LFBLE+LCITeJ)/CFLOATCIC)) INLET 589 

94 BCLs JJSRX(NI oT 9 J) INLET 390 
95 CONTINUE INLET 394 
CALL LEVEL INLET 392 
ASEN, INLET 393 
AHEN, INLET 394 
AF=0, INLET 395 

DO 97 TeleIs INLET 396 
AAZ0, INLET 397 
OLE0, INLET $9R 

00 96 JaleIC INLET 399 
DOLEML#L (Te JI/(CFLOAT(IC)*LE) INLET 400 

D(C LeJISNX(NIo le J) FHC I ed) INLET 40} 
IFCO(C Lod) ol 100.) OC Tod) 800001 INLET 402 

ACTe JI SACI op J)*OCI OJ) +HC Tod) #ABSCHC Lod) SC ZETAMFLOATCIC)) INLET 403 
IFCACT oJ) ob Tee) ACTS) #00001 INLET 4ou 
IF(T,EG,1) ASEAS+A(T oJ) INLET 405 


38 


% 


97 


159 
140 


64 
83 


85 


C FIND THE RELATIVE RANK MF TERMS, NORMALIZE BY THE LARGEST TERM, 


IF(T,EQ,1S) ABSABtA(T od) 

AAZEAA+A(T OJ) 

IF (AdoLToAMIN(NI)) AMINONI) BAA 

AEZAE+DL/AA 

AMINTONT)SAMIN(NT) 

AEB SAE 

IF(THT,EG.1) CALL WT} 

IFCIMT.EGo2) CALL wt2 

IFCYwT,FQ,5) CALL wt 

0O 140 Ysi,1S 

00 439 Jaie,Ic 

HX(NIoTeJ) BH(1eJ) 

wX(NToTeoJ) a¥*(TeJ) 

CONTINUE 

RNK(NIe 2) FAEs/(CO#LEV (Lo /(ARFR2) O1,/(AS#H2) ) *0Q*OO 
RNK(N]T 93) SG#AE/LE* (HBeHS) 

00 85 Talets 

ACEn, 

OO Aw JeloJC 

ACEAC#A(T oJ) 

DO AS JeleIC 

RNRONT eG) BRNKONTO4) SAF /CLE*AC)#GEN( Tod) #8 2M ABS (W(T oJ) #00)* 
PH (Te J) *ROS(2e20A8D(T 9 J) FF oS SSSSRA(T oJ) me2) aL (Tod) *B(IeJ) 
CONTINUE 

RNK(NI91)SeRNK(NJ 02) eRNK(NI gS) @RNK(NI 04) 
DERVONI)ERNK(NT01) 


XMAySO, 

DO 10! T#e1,4 

IF CABSC(RNK(NIOT)) oGT,XMAX) KMAXSABSCRNK(NI9T) ) 
ON 102 Tet,4 

RNK(NT eT) 8100 ,#RNK(NIo I) /XMAX 

CONTINUE 

DERV(NINLET#1) EQT/ARAVPOINFLO/SABAY 

RETURN 

END 


SUBROUTINE TRPwWRTECNINLEToXeHSoGINFLOOVo AMINT oO RNKoNT) 


C SUBROUTINF TO RITE HYDRAULIC INFORMATION ON TAPES 


DIMFNSIONN RNK(304) 9¥(5) oAMINI(3) 
HOURSEX/3H00. 

NTSNT¢) 

DO 100 NI®&{oNINLET 

IUNTTSNI¢A : 

VeYC(NI) ZAMINI CNT) 


WRIFECTUNIT) HOURS» HSeQINFLOVOYCNINLET#I) OVO YONI) 0 CRNK(NI od) oJB$ 04) 


RETURN 
END 
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INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
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INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
TNLET 
INLET 
INLET 
INLET 


406 
407 
408 
409 
410 
it 
412 
4\3 
414 
415 
416 
4\7 
418 
419 
420 
Wes 
422 
423 
4eu 
4es 
“eo 
Ue7 
428 
429 
430 
431 
432 
433 
454 
435 
456 
437 
438 
439 
440 
aus 


442 
4443 
4ua 
4us 
4U6 
uu7y 
44uB 
4ug 
490 
451 
492 


SUBROUTINE LEVEL 


C THIS ROUTINE COMPUTES WATER LEVELS THROUGHOUT THE INLET ASSUMING LEVEL 


C ARE LINFAR FROM BAY TO SEA 
REAL LeLENGTHeLINe LX eNoNX 
COMMUN/NUMS/NI gp GoNINLETOICH(3) eo ISE(3) eQGRokL (707) eB( 707) DC 797) 0 


1 80707) oNC707) oN 0797) 0V6707) 000797) eHSeHBoh( 797) eICeISeAMINI(3) 0 


{BMINTCS)eLINeQX(3) pQINFLOeABAY LENGTH(3) 
00 20 JeleIC 
XLEo, 
DO 10 Ieleys 
10 XLEXL*¥L(IeJ) 
XXEL (le) 726 
Hl Le J) BHS¢ (HBEHS) /XLEXX 
0O 11 I#eels 
KXT(LCLelLeJ)¢bL (195) )/2—9%XX 
11 HCIe J) 3HS* (HB&HS) XL OKX 
20 CONTINUE 
RETURN 
END 


SUBROUTINE SEACHSe TIME) 
C THIS SURRUUTINE DETERMINES THE FORCING SEA LEVEL EITHER FROM 
C EQUAL®TINE@SERIES DATA (IF AVAILABLE) OR BY SINUSODIAL FORCING, 
CUMMON /NUM3/A0g Te ARG BETA 
DIMENSION y(52) 
NNENNS] 
IFCNNeNE.1) GO TO 10 
REAN(S5e1) TOELeNPTS 
1 FORMAT (34X oF Oe206X_73) 
TOELETDFL¥*60, 
C READ SFA LEVEL EGUAL TIME SERIES DATA YHE FIRST TIME SEA IS CALLED 
C IF NPTS IS GREATER THAN 1 
IF(“pPTS.GT,3) READ(Se2) CYC J) oJB@leNPTS) 


2 FOKMAT(AF 10.5) 

IF(NPTS.GT,1) WRITEC693) (Y(J)oJalgNPTS) 
3 FORMAT( 3X01 6F 6,2) 

NISNPTSe#]1 

N2ENPTS¢2 


YONY)BYC1) 
y(neyey(e2) 
10 IF(NPTS,LT.1) GN TO 100 
C INTERPOL aTE IN TIME 
ITSTIMEST 
XTETINEmIT#T 
JEXT/IDEL 
JaJde4 
HSBY(J)o (Cy (JFL eV (J) eCXTe( Jel) *TOFLI/TUEL ) 
RETURN 
C DETERMINE LFVEL IF SEA LEVEL FLUCTUATION IS SINUSODIAL 
100 HSEAN* SIN(2,*3,14158"TIME/T) 
RETURN 
END 
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453 
45a 
455 
456 
457 
458 
459 
460 
“oi 
462 
463 
4ou 
465 
466 
4o7 
4o8 
469 
470 
471 


472 
473 
474 
475 
476 
477 
478 
479 
460 
461 
4b2 
483 
dou 
465 
466 
467 
468 
469 
490 
O91 
492 
493 
494 
495 
496 
497 
498 
499 
500 
501 
502 
503 


SUBROUTINE HELMCTHELMe AK CORL) 
C ESTINATE THE INLETe@RAY HELMHOLTZ PERIOD 
C OF THE TNLET/BAY SYSTEM (NEGLECT FRICTION) 
REAL LeLENGTHeLINOLX eNoNX 
COMMON/NUNG/NI gp GoNINLETOICH(3) 0 ISE(3) eOROL (707) eB( 707) oD 797) 0 
1 ACTe7) oC 707) ow (707) 0V6C797)9A(7T07) oHSoHBoH( ToT) oICoISHAMINI( 3) 0 
YRMINICS) oLING QX03) oA INFLOOARAY OLENGTHCS) 
OIMENSION CORL(3) 
C USE FIVE ITERATIUNS TO OBTAIN THE ESTIMATE 
00 1000 IIs195 
SUN=S0 
DO 100 NNZ1eNINLET 
AMINZAMINI (NN) 
100 SUHZSUM+AMIN/CLENGTHO(NN) ¢CORL(NN) ) 
THELME2,%3,1415A% SARTCAB/G)/ SART(SUM) 
C ESTIMATE THE HELMHOLTZ PERIOD 
DO 401 NNe{eNINLET 
C ESTIMATE THE INLET LENGTH CORRECTION DUE TO RADTATION 
101 CORLONN)S@BMINT CNN) 3, 141588AL0G(3,1415R*BMINT(NNI/S( SORTC 
13202*AMINI (NN) SAMINI (NN) )*THELM) ) 
1000 CONTINUE 
C CONVERT THE HELMHOLTZ PFRIOD TO HOURS 
THELMETHELM/36000 
RETURN 
END 


SUBROUTINE WTI 
THIS SURROUTINE WEIGHTS THE FLOW IN EACH SECTION SO THAT FRICTION 
IN THAT SECTION IS MENTMTZED, THIS MEANS THAT AT EACH SECTION FLOW IS 
ALLOWED TO REDISTRIABUTE ISTSELF THROUGHOUT THE CHANNELS TO MINIMIZE FR 
HOWEVER» FLOW PERPENDICULAR TO THE CHANNELS IS ASSUMED TO BE SMALL AND 
FLOW IS NOT INCLUDED IN THE EQUATIONS OF MOTION, BY MINIMIZING FRICTI 
ROUTINE GIVES aN UPPER LIMIT FUR RAY LEVEL FLUCTUATIONS AND INLET VELO 
REAL LeoLENGTHeLTNoLX oNo NX 
COMMON/NUMS/NI ep GoNINLETOICH(3) op ISE(3) eGAoL (797) eBC 707) 0D 797) 0 
1 AC7e7) ONC 797) 9H0797)0V 6797) 990797) oHSeHBoH( 797) 0TCoTSoAMINI(3) 0 
LAMINT CS) oLING QXC3) oQINFLOOARAY pLENGTH(3) 
DIMENSION C20) 
00 100 rai,Is 
SUMC=0, 
CO 50 Jeteo1C 
COJSACT oJ) *#2e (DCT) *%,3353)/ 
1 (NC T ed) **2"QX (NT) ewewB (led) #L(TeJ)) 
50 SUNCBZSUMCH+C (J) 
DO 60 JelelIC 
60 w( 19 J)8C(J) /SUMC 
100 CONTINUE 
RETURN 
END 
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504 
505 
506 
507 
5046 
509 
510 
Sid 
Sila) 
513 
Sig 
$15 
S16 
517 
518 
$19 
520 
Sel 
Se2 
523 
Seg 
525 
526 
S27 
528 


529 
530 
S51 
532 
553 
Sha 
535 
536 
537 
558 
539 
540 
Sui 
542 
543 
544 
545 
5u4 
547 
546 
549 
550 
551 


aaann 


aanaanananna 


SUBROUTINE we 


ROUTINE TO nETERMINE THE GRID WEIGHTING FUNCTION ASSUMING THAT 

FLOW IN A GIVEN CHANNEL IS THE SAME ALONG THE ENTIRE CHANNEL 

FLOW TS PISTRIBUTEN IN CHANNELS TO GIVE A MINIMUM TOTAL FRICTION 
FRICTION It) THIS ROUTINE WILL BE SLIGHTLY HIGHER THAN IN TL AND THE 
IN THIS SYSTEM ITS CONSISTANT WITH THE EQUATIONS OF MOTION, 


50 


60 
70 


1 


REAL LeLENGTHeLINeL Xo NoNX 


COMMUN/NUMS/NT 9 GoNINLETOICH( 5) oISE(3) oGAeL (707) 986797) 000797) 9 
1 ACTe7)eNC707) oW(797) 0V6707)90( 707) oHSeHBoH(Te7) 9 ICOISeAMINI(3) 0 
LHMINIC 3) oLINe OX (3) pM INFLOSABAY »LENGTH(3) 


DIMENSION C620) 
siimced, 

DG 100 TH1,9IC 
C(1j#0, 

00 50 JsileJsS 


COTyeC (1) + (NC Sol) Be 2eQX(NI) Fee (BC Jol ab (det))/ 


CAC Jel) # R24 ( 
CClys1,/C(1) 
sumceSUMC+C(J) 
DO 70 JeleIs 
00 40 TeleIC 
w(Je1)2C(1)/SUMC 
CONTINUE 
RETURN 
END 


SUBROUTINE wT3 


THIS ROUTINE ASSUMES THAT DISCHARGE 
THE INLET GRIN SYSTEM, JN GENERAL 
OIFFICULT TN ACCURATELY DRAW THTS TYPE OF GRID RY EYE AND FLOW DISTRUB 
CHANGES WITH TIME IN MOST INLETSe 
VELOCITTES AND BAY LEVEL FLUCTUATIONS, 
GRINS WITH NEPTRS LT 0,01 FOOT ARE 


100 


rw 


REAL Lol ENGTHOLTNe LX 9NoNX 
COMMONSNUMSSNIe GoNINLET®OICH(S) 9 


D(Jo1))**,33333)) 


IS EQUALLY DISTRIBUTED THROUGHOUT 
THIS WILL NOT BE TRUE BECAUSE IT IS 


THIS ROUTINE IS USEFUL IN GIVING AN 
ASSUMED TO HAVE NO FLOW 


ISEC3) oOROL (707) 0BC 707) 9O(T07)0 


1 ACTe7) eNO TTI eN( 797) oV6IT97) 9Q0707) oHSeHBeH( To 7) 9ICOTSeAMINI( 3) 9 
LEMINT (3) ol IN9 OX (3) eQINFLOPARAYOLENGTH(3) 


| 


00 2 IleieIs 

xzIc 

DO 4 Jatelc 

IF (M( Led) ob T0901) Kexele 
TF(X,LE.0,) WRITE( 60100) NteIS 


FORMAT(///¢5Xe( ERROR == INLET HAS ORIEN UP Ag INDICATED IN WTS C04 
5k, CINLET@ (eo T4ef SECTIONw[olI4oss/) 


IF(x,LE.00) STOP 

00 3 JeieIc 

WOl yd) a1 o/X 
IF(OC1lod) ol 700,01) w(Tod)B0, 
CONTINUE 

RETURN 

END 
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552 
553 
554 
555 
556 
557 
558 
559 
560 
Sot 
562 
563 
Sou 
565 
566 
567 
568 
569 
570 
S71 
S72 
573 
574 
575 
S76 


S77 
578 
579 
580 
561 
562 
56% 
Sbu 
585 
566 
587 
568 
589 
590 
591 
592 
593 
$94 
595 
596 
597 
596 
$99 
600 
601 


SUBROUTINE TABLE 


C ROUTINE T0 wRITE A TABLE nF INSTANTANEOUS HYDRAULICS 


59 


ol 


100 


REAL LoLENGTHeLTNoLXyNeNX 


COMMONSNUMNS/NI oe GoNINLFTOICH(3) oe ISE(3) oO ROL (707) BC 707) 0D(Te7)9 


L ACTeTI ONC 7TE7) eNO 707) 0V6797) 9A 707) oHSeHbeH( ToT) oICoTSeAMINI (3) 0 


LHMIN TOS) oL ING OK 63) >QINFLOOARAVYoLENGTH(3) 
CONMONJNUMY /Y¥ 05) eDERY (5) eXONToOITWT eo ZETAGHH 


COMMON/NUM2/BX(3e707) oOX(C 307 OT) OHK(BeTeTI OWN (SeToTIOLX( 30797) ONX(3S 


1e797) 

COMMONZNUMUSRNK (304) 

DIMENSION NAME (4) 

DATA NAME/GHV(F PS)  HHACFT2) ecOHWEIGHT voHLEVEL 
HRS=X/35600, 

WRITE (C601) HRS 

FORMAT (1X0 (pom eeetceeee 2 ees ee cesescaseeen (7/ 
1SXe(TIMF, HOURS 3[,F6,3) 

DO 100 NIZ{eoNINLET 

VRITE( 6010) NIeHSeHRe Y(NI) 


FORMATCs910Ko TINLET (oT Se/010X0 (SEA LEVEL OF Ta (oF 7oco/o10Xe (BAY LEVE 


LLeF Te leFTeQe/esOXe (NISCHARGE oCFSmaleE10 edo /o2X_ (CHANNEL 
5 


{10 4 2 3 4 
ICS TCH(NT) 

ISETSE (HI) 

OO wv Jeyolc 

Of % Tstets 


6) 


SECT 


ACT oe JIEAX(CnT elo J) *ONX(NI oT oS) FHX(NT go L0d) )tHUCNIOTOJ) #ABS(HK (NI oT od 


1) )/CZETA®FLOATCIC)) 

IFCACTo I) ob T.0.01) AC ToJ)E0, 
VCTeJISV(NT)#@HKX CNT o Te S)SACI oJ) 

IF CACTI oJ) LEo%.01) Vi Ted)=0, 

TFC JEW 1) WRITE( 6950) JoNAMECU) OCHX(NIo Tod) 9 T@iolS) 
wRITE (6069) 

FURMAT(/) 

WRITE (6050) JoNAME(1) o (VC Tos) oImledS) 
FORMAT(UX oT 20 5x eAbePKe OF 1002) 

WRITE(CO050) JoNAMF(P) oe (ACIOS) 0181 07S) 
*RITE(O050) JoNAME(3)o(MXCNTo9lod) ol@lols) 
CONTINUE 

WRITE (6959) CRNK(NTo IT) oTI@4 94) 


FORMAT(SXe (TEMP ACCE (oF Tole ( CONV ACCB(oF7,10( HEADBloF Tole l FRICS 


1(oF 7,1) 
VBARSY(NI) JAMINT(NT) 
WRITE (6061) VBAR eAMINI(NI) 


FORKMAT(5X0 (MEAN VELOCYTTY AT THE MINIMUM AREA SECTION@ (oF 7o20l FT/S 


1ECte lf AMINE (oF Oo2ef FTC() 
CONTINUE 

RETURN 

END 


43 


INLET 
INLET 
INLET 
INLET 
INLET 
IMLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 


602 
603 
604 
605 
606 
607 
608 
609 
610 
611 
612 
613 
614 
615 
616 
617 
618 
619 
620 
621 
622 
623 
624 
625 
626 
627 
624 
629 
630 
631 
632 
633 
634 
635 
636 
637 
638 
639 
640 
641 
642 
643 
64u 
645 
646 
O47 
648 
649 


aan 


2012 
2014 


SUBROUTINE CRITCNTeNELTOIUNTTeTeoNCYCLES) 


SUBROUTINE CRIT COMPARES 3 CONSECUTIVE FUNCTION POINTS 
AND WRITES MIDDLE POINT IF IT IS A CRITICAL POINT 


DIMENSION F(305) eMARK(5) oTERM(U) 

DATA MARKA/sSIH />e MARKB/IH®/ 

REWIND JUNTT 

NLINES=0 

TE=Ts/3600, 

eRITE (601009) 

00 3 Nzqe2 

REANCIUNIT) XeCPCNe J) eUFLeSVOCTERM(TI oT 18104) 
00 100 Ns3,NT 

REANCIUNIT) KeCF(3e J) oSBLOS)OCTERM(II) oTI1B4 04) 
IF(x LT.=1,.9E*1N) Gn TO 101 

youTs0 

OO pueoO 14 & le S 

PARK (TA) 3 MARKA 

IF (FCQeTA) © F(101A)) 20325 2020+ 2014 

IF (FC3eTA) = FCQeTA)) 20P0, 20150 2015 

IF (FC 3014) @ F(20TA)) 2015¢ 20150 2020 


C CRITICAL POINT VALUE FOLIND 


2015 


1111 
2101 


1009 
1 


1011 
1 
i 
1 


your s 

MARK(IA) @ MARKH 

IFC TAeEQ,1, AND, F(29TA) gGT,0,) HSHSF (201A) 
TFCT&eEN 1 gAND,F(2yTA) oGTL0,) T1EX 

TFC TASER G1 AND, F(OyTA) LEO.) HSLEF(2VTA) 
TFC TACEQD a1 gAND,F(CeTA) LEO.) TEx 
IFCTASZEQN, 3, AND, F(3eTA),GT.9,) HRH=F (Se A) 
TFC TACENs Sg4ND,F(SeTA) eGT20,) T38X 

IFC TAGENS SAND FC Se TA) LE GV,) HHLBF( SOTA) 
IFC TASER S,AND FC 3p TA) LE V,) TUBX 
IFCTASEN 4 AND F(OeTA),LT.O,) VEBF(20TA) 
IF (TAeED ed AND, F(2eTA) gGT00,) WEBF(QeTA)y 
CONTINUE 

UG 202d IA & te S 

F(CleTA) 3 Fl2eTA) 

F(2,TA) @ F(3eTA) 

IF (IGUT.EQ.0) GO TO 100 

IFC x,LT,(NCYCLES=2)"TF) GO TO 100 
NLINESZNLINES#{ 

IFCNLINFS.GT.150) GO TO 100 


WRITE (6 2101) Ko(FC{eTA)oMARK(TA) eo TAB105) 
CONTINUE 

RENT 

AMPHEHBH/HSH 


AMPL@HBI /HSL 

PHH=e AKS(T3"1T1)*360,/TF 

PHL= ABS(TUeT2)*360,/TF 

WRITE (691011) AMPHyPHHe VFO AMPL eo PHLeVE 

WRITE COoL111) TF 

FORMAT( 5Sxe(TFol9F7,2) 

RETURN 

FORMAT (2FBe3Se Alem SPFAeSVALVC(OPFTg 3HAl)9 

BPFQ,35 Aly 2CFA, 3S, A1)) 

FORMAT (UXeGHTIME + SX eQHHSeUXeOHINFLOW OSX e2HHRY 
Sve SHVEL se TX LHGe/eS5Xe SHHRS OSX GCHFT 9 SX gp U4HKCES 9 
OX ee2HFT 9 SX oe 3HFPSe HX oe UHKCESH/) 

FORMAT(C//91X9 (* CRITICAL POTNT VALUE [o///915X0 


(WAVE PROPAGATION (¢/915X0 (AG/A0 (05X09 (PHASE LAG(DEG) 


J/0 2K (HIGH WATER, 2K eo 3FI0eUesy 
2X_ (LOW WATER [(92xX93F10,4) 
END 
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650 
651 

652 
653 
654 
655 
656 
657 
658 
659 
660 
661 

662 
663 
664 
665 
666 
667 
668 
669 
670 
671 

672 
673 
674 
675 
676 
677 
678 
679 
6860 
661 

682 
683 
664 
645 
666 
687 
688 
669 
690 
oot 

692 
693 
694 
695 
696 
697 
698 
699 
700 
701 

702 
703 
704 
705 
706 
707 
708 
709 
710 
711 

Tle 


an 


10 
20 
25 


30 


50 


annan 


iselelelelelalelelelelelelelelelalealal 


SUBROUTINE READIN (Xe Vo YFAC XFACoX0oXF ep TNOCeKKeLNOIUNIT) 
SUBROUTINE TO READ SOLUTION TABULATION FROM FILE 


DIMENSION ¥(9)—9 YFAC(9) 
OT5S=,5*1./60, 

REAN CIUNTT) Xe Y 
IF(y,LT.91,E%10) KKae 

INDC ®@ 06 

IF (KK @ 1) 105 10, 50 

IF (X90 © X = DTS) 205 509 50 
IF (x © XF @ DTS) Joo 259 25 
KK go 2 

GO TO 50 

INDC ® 4 

X 8 XFAC#(X © XO) 

YCLN) B YFAC(LN)*Y(LN) 
RETURN 

END 


SUBROUTINE GRPHCC(ALABLIoALARL2>DELT»IUNTToNT) 


SURROUTINE GRPHC WRITES PLUTTER TAPE FOR GRAPHICAL 
OUTPUT OF SOLUTTON 


DIMENSION BL(2)> ISyHes) 

DIMENSTON YLARLL (3) oALEGNO $96) eALABLIS (4) oALABL2C4) op S¥M(3) 0¥(9) 
1669) oXX (2000) evy (20ND) eTT(9ee) 

DATA YLABLL/{OHHEIGHTS»s VolonHELOCITIES@yBHeFT» PRS/ 


oYPA 


DATA ALFGN/1OHFLOW ¢(KCFSe10H) 03H eLOHINLET VELO,1OHCITY 
1 (PT/Se3HEC) eo 10HBAY LEVEL (91 0HFT) 03H ol OHITNELOW o10H 
2 e3H eS OHOCEAN LEVEe1OHLC(FT) 03H oe lOHLEGEND o0H 
3 03H / 

DATA BKL/1OHOKSERVED Be lOHAY TIDE ( 


DATA ISYM/SeUodeevty/ 
DATA TT COL) /10HTEMPORAL A/ 


DATA TT(602)/10HCCEL / 
DATA TTC Tel) /10HCONVECTIVE/S 
OATA TT(702)/10H ACC / 
DATA TT(Be1)/10HPRESSURE H/ 
DATA TT(B92)/10HEAD / 
DATA TTC9e1)/10HBOTTOM STR/ 
DATA TT(902)/10HESS / 


READ INFORMATION TO DIRECT PLOTTING 


FIRST CaRD 
XQ = STARTING TIME OF PLOT (HRS) 

XF = ENDING TIME OF PLOT (HRS) 

SCALX = TIME aXIS SCALE TN HOURS PER INCH 

YLO = MINIMUM VALUE OF TIDAL HEIGHTS (FT) 

YL © OVERALL HEIGHT OF PLOT (INCHES) 

YLSC4L = SCALE OF TIDaL HEIGHTS (FT/INCH) 

YRO ® MINIMUM VALUE OF FLOWS (THOUSANDS OF CUBIC FEET PER SECOND) 
YRSCAL e SCALE OF FLOW ( THOUSANDS OF CUBIC FEET PER SECOND/INCH) 


CARD 2 

YVO © MINIMUM VELOCITY (FT/SEC) 

YVScalL = Scale OF VELOCITY (FEET PER SECOND/INCH) 
SCAL = SCALE FACTOR FOR TOTAL PLOT SIZE 

IQ e NoT EQUAL TO ZERQ FOR A PLOT OF INLET DISCHARGE 
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713 
714 
715 
716 
717 
718 
719 
720 
721 
Tee 
723 
724 
725 
726 
727 
728 
729 
730 


731 
732 
733 
734 
735 
756 
737 
738A 
739 
740 
74h 
742 
745 
744 
745 
7146 
747 
748 
749 
750 
751 
752 
753 
754 
795 
7196 
757 
758 
759 
760 
Toh 
7To2 
763 
Tou 
7605 
766 
7o7 
708 
709 
770 
aT 
772 


2001 


2002 


C DET 


1215 
C 
C 
Cc 


20 
C PLO 


IF (NITEM) 

TREAD (€ 502001) xOeKFe SCALXoVLOeYLeYLSCAL eYROVYRSCALeYVOOYVSCALY® 
1 SCaLEsIQ 

FORMAT (RF 10059/03F10.50110) 

KRITE( 602002) KOK eSCALX*YLOPYLeYLSCALPYROVYRSCALeYVUSYVSCALe 
1 SCALE+I0 

FORMAT(C///65X0 [PLOT INFORMATION[0/ 

1°1XeBF10eSe/4e1Xe3F10,50110) 

ERMINE SYMBOL SPACING 

LINTYP=,25*SCALK/(DELT/36000,) 

WRITE (CO01215) LINTYP 

FORMAT(1Xe (LINTYPS (916) 


PLOT LEGEND 


CALL SYMBOL (1 oeeYL/2,%7_80 02H OHLEGEND 00,06) 
DO 20 LN & tly S 

INDx = 06 

YPEovL/2,e, hel Nee2 

LLNSISYM(LN) 

CALL SYMBOL (Ose YP+ Oboe l4eLLNoOeenl) 

SYM(1) 3 ALEGN({9LN) 

SYi4¢2) 3 ALEGN( ALN) 

SY#(3)5 ALEGN(3eLN) 

CALL SYMBULC ed4eYPoMgt oSYMeVeeed) 

CONTINUE 

Nemlelan tee 

CALL SYMBOL C5e5e"YL 2,21 ee e2lvALABL1 90,932) 
CALL SYMBOL (6 5e5e2YL 72,21 049,21 ALABL200,932) 


C PLOT AxXFS 


1234 


60 
c IF 


YLOz=YL/2e*YLSCAL 

CALL AXTS(Oce*ViL/2ee1HHVELUCITYs FT/SECy 160 ¥L090,9YVO 
1eYVSCAL) 

CALL AXTS(eeBveYL/2,e11HHEIGHTSe FTeollevbe%eeYLOvYLSCAL) 
CALL AXTSCOge*YL/2,e9HTIME® HRS0@G9 (XFHKX0) /SCALKX 00,0009 SCALK) 
IFC TQ NEO) 


SCALL AXIS((XF=XN)/SCALKOPYL eye 1LOHFLUM) KCFSem1L00YL #9000 YL/2e¥%YR 


1SCALsYRSCAL) 
IFCTG6EN.9) CALL PLOTCCE XF#KX0) /SCALK9*YL72,03) 
IFC TQeE M20) CALL PLOT(CKFRXO)/SCALXOYL/2002) 
CALL PLOTC(XF"X0)/SCALKOYL/2903) 
CALL PLOT(OaeYL/2002) 

YFAC(1) 2 Y./YLSCAL 

YFAC(2) = 0.001 /YRSCAL 

YFAC(3) = YFAC(1) 

YFAC(4) = 1e/YVSCAL 

YFAC(S) = yYFAC(A) 

DO 4234 ITe609 

YFAC(I1)3,003 

KFAC B® 1,/5CALX 

DD AS 1% 19 9 
In=0 00 NUT PLOT DISCHARGE 

IP (TGLEN,0,AND,T,EQ,5) GO TN BS 
CORSYL/2e4( 195) #008 

CALL PLOT (%,0 Og 3) 

KK = 1 

ISUAaO 

REmINO TUNIT 
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INLET 


INLET 
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INLET 
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INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
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INLET 
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INLET 
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INLET 


773 
774 
775 
716 
777 
778 
719 
780 
781 
762 
783 
764 
785 
786 
787 
788 
789 
790 
uM 
792 
7193 
794 
795 
796 
797 
798 
7199 
660 
Bol 
802 
603 
604 
805 
806 
507 
608 
609 
610 
811 
812 
613 
614 
615 
816 
817 
818 
619 
820 
B21 
Bee 
823 
824 
625 
626 
627 
828 
629 
650 


65 


70 
Te 


60 


C PLOT CURVES ( nO NOT PLAT IF EWUAL TO ZERO THROUGHOUT) 
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C READ PROTOTYPE BAY TIDE 


2019 


INDY 
CALI 


GO TU 


e494 


READIN (XeVeYFACeXFACoKOeXFoINDOC oKKe lo IUNIT) 


(7Ce BO)» KK 


TFC rNOC,LE,0) GON TO 65 


ISUAS 


TStibey 


IF (TSUb,GE,3998) 1511B21998 
XX (7TSUnyex 
YY( TSI) zYC(T) 


TEGyy 


GT,5) YYCISUB)syyY(ISUB)+COR 


IF(TSUB,EW,1998) GO TO 80 
GO TO 65 

XXCTSUKe1) 50, 
MX(TSUbe2)51,0 
YY(TSUBK41)50, 
YYCTSUbe2) 21, 


TFC vy CISUBe2) eENe9 gM gANDo 


1 YY¢TSUR#@1),EQ,0,9, ANN,YYCISUA) oEQ,0,0) GO TO 65 
IFCY, 


CALL 


61,5) 60 TN BBS 
LINE CXXeVVYeISUReloLINTYPel) 


GO TO 85 


CALL 
CALL 
CALL 


LINE CKXeVVeTStiAg19090) 
PLOTC(XFex0)/SCALXeCOR, 5) 
PLOT(0eeCORe2) 


SYMCLISETTCI01) 
SYN(2)ETT(1Te2) 


CALL 


SYMBOL Ce2,20CORVDOeloSYMo0e020) 


CONTINUF 


TF CNITeNF,1) GO TO 2019 
REAN(Se1) TOELeNPTS 
FORMATC SUX» FoeeobXoT3) 
IF(NPTS,LT,2) GN Tn 2019 


IF(NPTS,GT,1) READ(So?) (CYY(J)oJ@1oNPTS ) 


FORMAT(AF 1005) 
XX(NPTS#1)2£9, 
XKX(NPTS#2)e1, 
YY(NPTS+1)50, 
YY(NPTS#2) 21, 


po 3 


Jay{oNpTS 


YY(CJ)SYV (J) *YF ACL) 
XXCIIECTOEL S69, )*XF ACH (Jel) 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 


PLOTCXX(1)0YY(1) 03) 

LINE CKXe YY eNPTSo10000) 
PLOT(CXX(NPTS/2) oYY(NPTS/2) 03) 
PLATCKXCNPTS/2) eYYCNPTS/2)4,7502) 


SYMBOL CXX(NPTS/2) +01 9YY(NPTS/2)¢e75001obLo00017) 


PLOT(C(XF@XO)/SCALK+4 400,095) 


RETURN 


END 
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(MATA STARTS AT BEGINNING OF PLOT» SAME NATUM) 
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INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
INLET 
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INLET 
INLET 
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832 
633 
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658 
859 
640 
B41 
842 
Buy 
644 
845 
846 
647 
648 
649 
850 
651 
852 
653 
854 
655 
856 
657 
658 
859 
860 
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B62 
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665 
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868 
6609 
670 
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872 
673 
874 
875 
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